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Abstract 

We present the up to now most precise evaluation of electroweak and super- 
symmetric contributions to the anomalous magnetic moment of the muon a^, 
describing in detail also the calculational techniques. We calculate the bosonic 
two-loop contributions in the Standard Model without the approximation of 
a heavy Higgs-boson mass, finding corrections up to 0.2 x 1CP 10 for a light 
Higgs boson. In the Minimal Supersymmetric Standard Model the correspond- 
ing two-loop contributions from the two-Higgs-doublet model part differ from 
the Standard Model result by up to 0.3 x 10~ 10 . Finally, we evaluate the dia- 
grams where a loop of charginos or neutralinos, the superpartners of gauge and 
Higgs bosons, is inserted into a two-Higgs-doublet one-loop diagram. These 
corrections can amount up to 10 x 10 -10 , which is almost 2a of the current 
experimental uncertainty. 
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1 Introduction 



The final result of the Brookhaven "Muon g — 2 Experiment" (E821) for the anomalous 
magnetic moment of the muon, a^ = (g — 2) M /2, reads [1] 

a e ^ = (11 659 208 ±6) x KT 10 . (1) 

The Standard Model (SM) prediction depends on the evaluation of the hadronic vac- 
uum polarization and light-by- light contributions. The former have been evaluated 
by Refs. [2-5], the latter by Ref. [6], but there is a recent reevaluation [7], describing 
a possible shift of the central value by 5.6 x 10~ 10 . Depending on which hadronic 
evaluation is chosen, the difference between experiment and the SM prediction lies 
between the two values 1 

a c xp _ a thco ^ pj + ^ = ( 3L7± 9.5) x 10 -io . 3 3 a ; ( 2 ) 

a ex P _ a th«, ([ 2 ]+ [7]) = (20.2 ± 9.0) x 1(T 10 : 2.1a . (3) 

These evaluations are all e + e~ data driven. Recent analyses concerning r data indi- 
cate that uncertainties due to isospin breaking effects may have been underestimated 
earlier [4]. For the purpose of numerical comparisons we will use the intermediate 
value 

q^p - a^ 60 ( [3]+ [7]) = (24.5 ± 9.0) x lO" 10 : 2.7 a . (4) 



It is an interesting question whether the observed 2 — 3a deviation is due to su- 
persymmetric effects. The supersymmetric one-loop contribution [9] is approximately 
given by 

a susY,iL = 13 x 10 - 10 /lOOGeVy ^ 

V M SUSY / 

if all supersymmetric particles (the relevant ones are the smuon, sneutralino, chargino 
and neutralino) have a common mass Msusy- Obviously, supersymmetric effects can 
easily account for a (20 . . . 30) x 10~ 10 deviation, if /i is positive and Msusy lies roughly 
between 100 GeV (for small tan/?) and 600 GeV (for large tan/3). This mass range is 
both allowed by present search limits and very interesting in view of physics at Run II 
of the Tevatron, the LHC and a future Linear Collider (LC). 

Eq. (jHJ) also shows that for certain parameter choices the supersymmetric contri- 
butions could have values of either a^ USY > 60 x 10~ 10 or a^ USY < — 10 x 10~ 10 , both 
outside the 3a band of the allowed range according to (0), ©• This means that the 
(g — 2)^ measurement places strong bounds on the supersymmetric parameter space. 
This is important for constraining different variants of supersymmetric models (and of 
course also other models of new physics) and complements the direct searches. Even 
after the discovery of supersymmetric particles, indirect bounds derived from (g — 2) /t 

1 We always include the updated QED result from Ref. [8]. 
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or 6-decays and a Higgs mass measurement will provide important complementary 
information to that obtained from direct measurements. 

In order to fully exploit the precision of the (g — 2) M experiment within Supersym- 
metry (SUSY), see Refs. [10,11] for possible applications, the theoretical uncertainty 
of the SUSY loop contributions from unknown higher-order corrections should be 
significantly lower than the experimental error given in eq. and the hadronic un- 
certainties in the SM prediction. Thus, the reduction of the uncertainty of the SUSY 
loop contributions down to the level of about ±1 x 10~ 10 is desirable. 

For the electroweak part of the SM prediction this accuracy has been reached with 
the computation of the complete two-loop result [12,13]. However, this result relies 
on a single evaluation of the bosonic two-loop contributions, performed in the limit 
of a heavy SM Higgs-boson mass, M h sm 3> My/- Thus, an independent check seems 
desirable. In this paper we perform this evaluation for arbitrary Higgs-boson masses. 

For the SUSY contributions, a similar level of accuracy has not been reached yet, 
since the corresponding two-loop corrections are largely unknown. Only two parts 
of the two-loop contribution have been evaluated up to now. The first part are the 
leading log (m At /MsusY)-terms of supersymmetric one-loop diagrams with a photon 
in the second loop. They amount to about —8% of the supersymmetric one-loop 
contribution (for a SUSY mass scale of M SU sy = 500 GeV) [14]. 

The second known part are the diagrams with a closed loop of SM fermions or 
scalar fermions calculated in Ref. [15], extending previous results of Refs. [16,17]. It 
has been shown in Ref. [15] that the numerical effect of these contributions could in 
principle be as large as 20 x 10~ 10 for suitable parameter choices. But experimental 
constraints on the lightest Higgs-boson mass [18-21], electroweak precision observ- 
ables [22-24] and 6-decays [25, 26] restrict the allowed parameter space. Thus the 
values of these diagrams amount up to about 5 x 10 -10 , except in rather restricted pa- 
rameter regions with non-universal sfermion mass parameters involving very disparate 
mass scales. 

In general, every diagram in the Minimal Supersymmetric Standard Model (MSSM) 
must contain one continuous line carrying the /i-lepton number. Thus, the MSSM di- 
agrams can be divided into two classes: (1) diagrams containing a one-loop diagram 
involving supersymmetric particles (with a jx or z> M line) to which a second loop is 
attached, and (2) diagrams where a second loop is attached to a two-Higgs-doublet 
model one-loop diagram (with a /i or ^ line). The QED-logarithms from Ref. [14] 
belong to the first class, and the fermion/sfermion two-loop diagrams from Ref. [15] 
to the second. 

The diagrams of the second class are particularly interesting since they can depend 
on other parameters than the supersymmetric one-loop diagrams and can therefore 
change the qualitative behaviour of the supersymmetric contribution to a M . In par- 
ticular, they could even be large if the one-loop contribution is suppressed, e.g. due 
to heavy smuons and sneutrinos. 
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In this paper we complete the calculation of the diagrams of this class. As the 
first step we calculate all pure two-Higgs-doublet model diagrams. This calculation 
is analogous to the one of the SM bosonic two-loop corrections (we carry out the 
calculation of the SM contributions as a special case), but it involves the additional 
Higgs bosons (but no SUSY particles). The second step comprises the diagrams with 
a closed chargino/neutralino loop. Their structure is similar to the one of diagrams 
with a closed SM fermion or sfermion loop, but they depend on a completely different 
set of parameters and show a complementary behavior. 

The objective of this paper is to describe in detail the calculational steps, which 
have already been used in Ref. [15], and to analyze the numerical impact of the newly 
derived SUSY contributions. The SM result is compared to the existing calculation. 
As for the results of Ref. [15], the new results will be implemented into the Fortran 
code FeynHiggs [27]. 

The outline is as follows. In Sect. 2 we describe the calculation, in particular the 
used regularization, large mass expansion, reduction of two-loop integrals, and renor- 
malization. For the two-loop QED corrections involving SUSY particles analytical 
results are presented. Section 3 is devoted to the SM contributions. The numerical re- 
sults for the two-Higgs-doublet contribution in the MSSM and the chargino/neutralino 
contribution are discussed in Sects. 4 and 5, respectively. Our conclusions are given 
in Sect. 6. 



2 Calculation 



2.1 Extraction of a. 



The anomalous magnetic moment a p of the muon is related to the photon-muon 
vertex function T^ap as follows: 

u(p')T^ AP (p,-p',q)u(p) = u(p')[ lp F v (q 2 ) + (p + p') p F M (q 2 ) + ...]u(p), (6) 



a p = -2m M F M (0). 



(7) 



It can be extracted from the regularized vertex function using the projector [12, 13] 



1 



D-2 



2{D - 1){D - 2)ml 1 2 



Tr\ — — [mj 7p - Dp P i -(D- l)m pPp ] V 



m, 



V p = T pflAP (p,-p,0), 
d 



(8) 

(9) 
(10) 



q=0 
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Here the muon momentum is on-shell, p 2 = m 2 , and D is the dimension of space-time. 

This projector requires that in the covariant decomposition (jUJ) only D-dimensional 
quantities appear. We therefore use dimensional regularization with anti-commuting 
75, where there is no distinction between the first four and the remaining D — 4 dimen- 
sions. However, in order to demonstrate the validity of this regularization we need 
to discuss three issues related to supersymmetry breaking, mathematical consistency 
and the treatment of the e-tensor in this scheme. 

• Dimensional regularization breaks supersymmetry [28], and in general super- 
symmetry has to be restored by adding certain counterterms not correspond- 
ing to multiplicative renormalization of the fields and parameters [29]. In our 
case, however, all appearing counterterms are two-Higgs-doublet model coun- 
terterms, because we calculate 2-loop corrections to 1-loop two-Higgs-doublet 
model diagrams, see also Sect. 12.31 These two-Higgs-doublet model countert- 
erms are either fixed by renormalization conditions, like the muon- or Z-mass 
counterterm, or by gauge invariance, like the /i/ry- or H /+ G~7-counterterms. 
Since gauge invariance is not broken by dimensional regularization with anti- 
commuting 75, 2 multiplicative renormalization is sufficient for all counterterms 
in our calculation. 

• Using an anti-commuting 75 in D 7^ 4 implies Tr(757 M 7* / 7 p 7 <J ) = and is 
therefore incompatible with the trace formula Tr(7 5 7 M 7 iy 7 p 7 <J ) oc e M ^ pfT , which 
has to be reproduced for D — > 4. As it is customary, we simply replace each 
trace of the form Tr(7 5 7 M 7 I/ 7 p 7°") by its four-dimensional value. We can check 
the correctness of this procedure by comparing the terms involving e-tensors with 
the ones obtained using the HVBM-scheme [31], which is fully consistent. In our 
calculation, e-tensors appear only in two-loop diagrams with an insertion of a 
fermion triangle with three external vector bosons. The difference of the fermion 
triangle subdiagrams in the two schemes (summed over one generation or all 
charginos and neutralinos) is of the order (D — 4)|f ^ v P a ; where k a is the non- 
vanishing external momentum of the fermion triangles. The covariant e^ p<T k a of 
power-counting degree +1 has the same prefactor as the chiral gauge anomaly 
and thus vanishes. If the fermion triangles are inserted into 1-loop diagrams 
contributing to T^^ap, the difference between the naive and the HVBM-scheme 
vanishes for D — > 4 for power-counting reasons. Hence, the terms involving e- 
tensors are identical in the naive and the HVBM-scheme, which means that the 
naive scheme produces these terms correctly (see also the discussion in Ref. [30]). 

• The e-tensor is a purely four-dimensional object, and contractions like e tia ^e v a p 1 
or e^V^ have to be evaluated in four dimensions. 3 Effectively, this leads to 
the appearance of the additional covariant y^-dim^i-dim - m ^ This con _ 
tradicts the requirement of the projection (jHJ) that only .D-dimensional covari- 

2 This has been checked up to the order we need here in Ref. [30]. 

3 We checked explicitly that a Z?-dimensional treatment indeed leads to incorrect results. 
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ants may appear in T^ap- However, as mentioned above, the fermion triangle 
subdiagrams add up to an expression of reduced power-counting degree, and 
accordingly the e-tensor contributions to T^ap and thus Fy are finite. As a 
result, the projection operator (JHJ) produces the correct result, containing only 
the form factor Fm, even if the e-tensors are evaluated in four dimensions. 

2.2 Diagram evaluation 

The actual calculation is done using computer algebra. The Feynman diagrams are 
generated using FeynArts [32,33]. After applying the projector (jHJ), the Dirac algebra, 
traces and contractions are performed by TwoCalc [34]. 

The main part of the two-loop calculation consists of the evaluation of the two- 
loop integrals and the simplification of their coefficients, both of which is complicated 
by the large number of different mass scales and the involved structure of the MSSM 
Feynman rules. 

As a first step we perform a large mass expansion [35] in the ratio m M /M heavy , 
where M heavy stands for all heavy masses, M Z) w an d Higgs and chargino/neutralino 
masses. Depending on the prefactors we expand all integrals to sufficiently high order 
such that the final result is correct up to 0(m^/M^ eavy ). We discard all terms of 
(9(m^/M^ eavy ). The expansion is based on the formula [35] 

Fr~^fr/ 7 or 7 F 7 . (11) 

7 

Here Fr denotes the integral corresponding to a (two-loop) Feynman diagram T, r/7 
denotes the diagram where the subdiagram 7 is shrunk to a point and 7^ denotes 
Taylor expansion with respect to all small masses and all external momenta of 7. The 
sum runs over all (in general unconnected) subdiagrams 7 that 

• contain all lines with heavy masses, 

• are one-particle irreducible with respect to the light lines. 

By means of (|llj) the large masses appear only in 7, whereas the small masses and 
momenta appear only in r/7, resulting in a separation of scales. Since the projection 
operator (jHJ) sets the external photon momentum to zero, our integrals are initially 
two-loop two-point integrals. There are therefore three possibilities for the r.h.s. of 
the large mass expansion (jll|) and the subsequent reduction to master integrals: 

• (Light 0-loop) o (heavy 2-loop): In this case, the r.h.s. of (fTT|) results in a two- 
loop vacuum diagram times a rational function of m^. Using partial integration 
identities [36], we reduce all two- loop vacuum diagrams to the master integral 

T 13A (m 1 ,m 3 ,m 4 ) = (( * — )) (12) 



5 



in the notation of [34]: 

A = " ml (13) 

h = qi, k 2 = qi+p, h = q 2 - q x , (14) 

h = q 2 , h = q 2 +p, (15) 

«■■■» = Swi^k^ (16) 

The result for the general case of three different masses can be found e.g. in 
Ref. [37]. 

(Light 1-loop) o (heavy 1-loop): In this case the integrals can be reduced to the 
standard one-loop functions A (m) and S o (m^,0, m M ) [38]. 

(Light 2-loop) o (heavy 0-loop): This case appears e.g. in the calculation of dia- 
grams 7,8 in Fig. below. In general it can only appear in diagrams that do not 
involve any SUSY particles, since SUSY particles (in scenarios with R-parity 
conservation) necessarily form at least one closed heavy loop. The r.h.s. of (jllj) 
then contains two-loop two-point functions, however with only one mass scale, 
m M . One typical example of such integrals is 

YiiH(ml;m„m„m,) = «®^». (17) 

jj 2 jj 3 u^ 

with p 2 = m 2 ^ and 777,2,3,4 = m M . Using Passarino-Veltman decomposition for the 
one-loop subdiagrams [34] and partial integration identities all such integrals 
can be reduced to one-loop integrals and the on-shell (i.e. p 2 = mV) master two- 
loop sunset integrals T 2 3 4 (m^; m^, m^, m M ) and T 2 u(Tn 2 l ] m M , 0, 0). For example, 
for (|17|) we obtain 

Y 2^!i m l^ "fy> "fy> "fy) = 

4(-4 + 19D - 31D 2 + UD 3 ) 4 2 

— = m,,A (m u Y (18) 

8-18D + 9D 2 K 1 

32(-2 + D)D 2 6 2 
3(8 - 18£> + 9D 2 ) ^ 

For the analytical results of on-shell two-loop two-point functions, see Ref. [39] 
and references therein. Note that the reduction can generate spurious additional 
divergences, in the form of 1/(D — 4)-poles in prefactors of master integrals. A 
simple example is given by the reduction of Y 2U5 with p 2 = m 2 and m 2 , 3j 4 = m M , 
m 5 = 0: 

*2345("4 «fy> "fy> "fy> °) = 

AA^m^) S o (mJ,0,m M ) 



"Hi 

1 



2(D-4) 

+ (£>-: 



^(wy) 2 
K(D-3) 



2T 2 3 4 (mJ;m M ,m M ,m At ) 



(19) 
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Although the 1/(D — 4) 3 -poles cancel, the analytical results of the master inte- 
grals are needed up to 0(D — 4) in the case of T 234 and up to 0((D — A) 2 ) in 
the case of A and B functions. 

2.3 Counterterm contributions 

Since a M is exactly zero at tree-level, no two-loop counterterm corrections arise. How- 
ever, one-loop diagrams with counterterm insertion are necessary to derive a UV-finite 
result. Thus, together with the two-loop diagrams we have also evaluated all the di- 
agrams with subloop renormalization in the SM and the MSSM. 

In order to derive the amplitudes of the diagrams with counterterm insertion it 
was necessary to extend the existing FeynArts MSSM model file [33] by including non- 
SM counterterm vertices. The introduced counterterms contain field renormalization 
constants for the gauge and Higgs bosons and for the muon and the neutrino. We 
have checked that in the sum of all amplitudes all field renormalization constants, 
except of the external muon and photon, drop out as required. In addition to the 
field renormalization constants also mass counterterms for the SM gauge bosons and 
the muon have been introduced. Finally, counterterms for the Higgs sector arising 
from the Higgs potential (see e.g. Ref. [20]) as well as for the mixing of Higgs- and 
gauge bosons have been included. The Higgs sector counterterms consist, besides field 
renormalization constants, of Higgs tadpoles that cancel large corrections coming from 
the self-energy insertions at the two-loop level. 

We have chosen on-shell renormalization conditions for the muon and the photon, 
for the W- and Z-boson masses, and for the electric charge. The tadpole counterterms 
cancel their corresponding loop contributions. 

Since all amplitudes are reduced to self-energy like corrections, the result of the 
counterterm diagrams consists only of A and B functions. As for the genuine two- 
loop diagrams, also for the counterterm diagrams a large mass expansion is necessary 
in order to obtain a consistent expansion in powers of m M . 

As already outlined in Ref. [15] we parametrize our one-loop result with G M in order 
to absorb process-independent higher-order corrections. This requires the evaluation 
of the one- loop corrections to muon decay, Ar, from the corresponding sectors. For 
the class of corrections evaluated in Ref. [15], all contributions to Ar from scalar 
fermion diagrams (since the SM fermions contribute only to the SM result) had to be 
taken into account. For the results computed in the present paper the Ar corrections 
from charginos and neutralinos, as well as from the whole two-Higgs-doublet sector of 
the MSSM (except fermions) had to be evaluated. While the former ones consist only 
of self-energy corrections, the latter ones also contain vertex and box contributions 
involving SM particles. The Ar contribution enters via (— a^ L Ar), where a^ L is the 
two-Higgs-doublet model one-loop result. In the case of the SM gauge and Higgs- 
boson contributions the Ar term amounts to a shift in the two-loop corrections of 
about 10%. 
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2.4 QED contributions 



In general, the contributions to a M can be split into QED and electroweak contribu- 
tions. A diagram is counted as a QED diagram if it involves only the photon but no 
other gauge or Higgs boson, and as an electroweak diagram if it contains a W, Z or 
Higgs boson (or some of the corresponding ghost fields). 

At the two-loop level, there are QED diagrams involving a photon vacuum po- 
larization subdiagram (see Fig. HJ. For leptons or quarks, these contributions are of 
course known; in the latter case they have to be obtained via the hadronic vacuum 
polarization [2-5]. 




Figure 1: Generic QED diagrams with a fermion loop insertion (left) and the corre- 
sponding counterterm diagram (right). The same kind of diagram exists with a scalar 
inserted instead of a fermion. 



If supersymmetric particles contribute, the two-loop QED contributions are mod- 
ified by the diagrams of Fig. ^ with a chargino, slepton, or squark loop. In our 
calculations in the following sections, as well as in Ref. [15], we do not include these 
QED contributions. They can be easily inferred from known results [40] and are tiny. 
For completeness, we list the results here. The result of these diagrams depends only 
on the charge Q, the mass M, and the spin of the particle in the vacuum polarization. 
For a scalar particle we have 

OED,2L,scalar _ Q 2 ( <*\ 2 K 



UEU,2L,scalar _ _^ / " \ M /r> n \ 

a » ~ 360 \7r/ M 2 ' 1 Uj 



and for a fermionic particle we have 



OED,2L,fcrmionic _ Q 2 ( a \ 2 m /i 



yED,2L,fcrrmomc ~e / \ /i /n-i \ 

~ 45 Vtt) M 2 [ZL) 

For masses M > 100 GeV, these contributions are below 10~ 13 and hence negligible. 

The results in f!20|) . (|2ip contain the suppression factor m^/M 2 , which is typical 
for all diagrams except for the pure QED-ones involving only photons and muons. 
However, this only holds after renormalization, i.e. for the sum of the two-loop dia- 
grams and the counterterm insertions. The two-loop diagram in Fig. ^ itself has a 
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contribution of the order m° which is then cancelled by the counterterm diagram. 
The reason is that before renormalization, the vacuum polarization subdiagram in- 
cluding the two photon propagators behaves like ■^IT' (k 2 ) oc for k 2 — > 0, like 
the tree-level photon propagator, whereas after renormalization (here charge and field 
renormalization effectively amount to adding the counterterm diagram in Fig. ^) we 
have ±rW(k 2 )± oc 1 for k 2 — > 0. 



3 Standard Model contributions 

Before discussing the MSSM results, we present our results for the bosonic electroweak 
two-loop contributions in the SM. Up to now, there exists only one evaluation of these 
diagrams [12], which employs the approximation M h sm 3> Mw- Our recalculation of 
these contributions serves both as a cross check of Ref. [12] and of our algebraic codes. 
It furthermore allows to compare the approximation of Ref. [12] with the result for 
arbitrary M h sm. 

Some typical diagrams of this class are shown in Fig. El below, if one identifies 
with the SM Higgs boson and ip with the charged Goldstone boson of the SM (in the 
Feynman gauge). In general this class consists of all SM two-loop diagrams without a 
closed fermion loop and without pure QED diagrams. We obtain the following result: 

HOS.2L = 5G^«/ ^ \ 

3 8ir 2 y/2 7T V Mfa J K ! 

where the coefficient of the large logarithm, log jtt- ~ —13, can be written in analyt- 
ical form: 

1 



c h L ° s ' 2L = ^ [107 + 23(1 -40 2 ] «3.6, (23) 



in agreement with Ref. [12]. The coefficient Co° s,2L has a more involved analytical 
form. In Fig. |2Jwe show the results for a^ os,2L from eq. (J22)) as a function of M h sm. 

The variation comes from the non-logarithmic piece Cq OS ' 2L . The size of the (constant) 
logarithmic contributions is also indicated. It can be seen that the variation with 
M h sm is at the level of 0.3 x 10~ 10 for 100 GeV < M h sm < 500 GeV. These numerical 
values also confirm the results of Ref. [12]. The logarithmic piece in (|22|) alone is 
already a very good approximation and leads to 

a bos,2L = (_ 22 ± Q ^ x 10 -io ; ( 24 ) 
corresponding to a reduction of the electroweak one-loop contribution by about 11%. 
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Figure 2: a^ os,2L is shown as a function of M h sm, see eq. 
stant) logarithmic contributions is also indicated. 



The size of the (con- 



4 Two-Higgs-doublet contributions 



In the MSSM, the bosonic electroweak two-loop contributions are different compared 
to the SM because of the extended MSSM Higgs sector. In this section we present 
the result for the contributions of this class, defined by selecting all MSSM two-loop 
diagrams without a closed loop of fermions or sfermions and without pure QED- 
diagrams, see Fig. El 



The result > os ' 2L ' MSSM reads 
a* 



bos,2L,MSSM 
V 



5 G^m 2 a 



3 8vr 2 V2 vr 



- c 



bos,2L,MSSM 



log 



M 2 W 



fj. . bos,2L,MSSM 

+ c 



(25) 



where the coefficient of the logarithm is given by 



bos,2L,MSSM 



C, 



^[98 + 9^ + 23(l-4 S 2 v ) 2 ] 



c 2/3 M 2 



111 



H 



m 



h J 



(26) 
(27) 



Here (3 is defined by the ratio of the two Higgs- vacuum expectation values, tan (3 = 
Vijv\\ nih^H and a are the masses and the mixing angle in the CP-even Higgs sector, 
and c a = cos a, etc. As we will discuss below, c£ = 1, and thus the logarithms in the 
SM and the MSSM are identical. 
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Figure 3: Some two-Higgs-doublet model diagrams for a M with (depending on the 
diagram) = h, H, A,G; $ = G ± , #±; V = <y,Z. 



Let us first comment on the non- logarithmic contribution in (J25j) . The non- 
logarithmic coefficient c Jj os > 2L ' MSSM j s a complicated function of the full MSSM Higgs 
sector, which is determined by two parameters, tan (3 and the pseudo-scalar Higgs- 
boson mass Ma- However, in Fig. 0] it is demonstrated that the full results for 

^bos,2L,MSSM q^-q 

be quite well approximated by the corresponding SM quantity if 
M h sm is set equal to the light CP-even Higgs-boson mass of the MSSM, M h sm = mh, 
except in parameter scenarios with extremely light CP-odd Higgs-boson mass. For 
Ma ^ 100 GeV the difference stays below 0.3 x 10~ 10 . The difference between the 
MSSM and the SM comes solely from the different coefficients c jj os > 2L ' M M and Cq OS,2L . 

As in the SM, the logarithmic piece in (J2"5j) is an excellent approximation of the 
full bosonic result. At first sight, the coefficient ()26|) seems different from the corre- 
sponding SM one, owing to the appearance of Higgs-mass dependent terms. However, 
the MSSM Higgs-masses satisfy the relation 



C /3 



CaC a +/3 S a S a -\-j3 
2 ' 2 
m H m h 



(28) 
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Figure 4: The difference of ajj os ' in the MSSM and the SM is shown as a function of 
Ma- The SM Higgs-boson mass has been set equal to the light CP-even MSSM Higgs 
mass, tan/? has been set to tan/5 = 2,50. The differences between the MSSM and 
the SM come solely from the different coefficients c b° s ' 2L > M ssM and Cg OS,2L . 



at tree-level. Hence, the coefficient & L = 1 in (|26|). and the logarithmic pieces of the 
MSSM and SM bosonic two-loop corrections are identical as mentioned above. This 
finding can be explained in several different ways. One way is to apply the analysis 

of the log TT7-- terms performed in Ref. [14] to the bosonic contributions. This shows 
that they must be identical in the SM and MSSM since the corresponding bosonic 
one-loop diagrams are identical. The latter statement is true since the additional 
MSSM Higgs bosons do not appear at the one-loop level, in particular since there is 
no 7iy ± if T -coupling. 4 

Alternatively, one can directly investigate the first Feynman diagram in Fig. |3] 
(with <p = h,H and ip = in order to understand how the logarithms in the SM 
and MSSM emerge. The diagrams with photon and physical Higgs exchange and an 
inner G ± — W loop are responsible for the Higgs-mass dependence in the coefficients 
of log jp-. In the corresponding SM diagram (with <\> = H SM ) the G + G~i/ SM -coupling 

is proportional to M^- sm /Mvk, and the dependence on M h sm exactly cancels with the 
suppression coming from the Higgs propagator. In the MSSM the couplings of G to 
physical Higgs bosons are given by gauge couplings, hence leading to the structure of 

4 More precisely, the additional MSSM Higgs bosons as well as the SM Higgs boson do appear in 
one-loop diagrams, but their contributions are suppressed by two additional muon Yukawa couplings 
and hence negligible. 
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c*l shown above in eq. ()27|) and to the seeming ^-dependence of the logarithm. 

However, as the appearance of Goldstone bosons signals, such Feynman diagrams 
are not gauge independent by themselves. For example, in the unitary gauge there 
are no Goldstone bosons and such diagrams do not exist. Similarly, in the background 
field gauge [41] or in the nonlinear gauge used in Ref. [12] there is no 7lV ± G =F -vertex, 
and again the diagrams in Fig.Eldo not exist. In all these gauges there are no diagrams 
with virtual photon and physical Higgs, hence it is manifest that the logarithms in 
the SM and MSSM are identical. 

In the Feynman gauge, which we we have chosen, this fact is not immediately 
obvious but follows from relation ()28|) , which in turn is a consequence of global gauge 
invariance. In fact, one can easily derive that the Higgs mass matrix T^.^ (corre- 
sponding to the denominators in (J28|) ). where faj are the CP-even interaction Higgs 
eigenstates, is related to the couplings Y^o+g- (corresponding to the numerators 
in (J2E1)). 5 

This discussion has a practical implication concerning the treatment of higher- 
order contributions to the Higgs-boson masses. In other sectors like the sfermion- 
loop or chargino/neutralino-loop sector, leading three- loop effects can be taken into 
account by using the loop-corrected values for the Higgs-boson masses M^h in the 
propagators. In the present sector, however, using loop-corrected Higgs-boson masses 
would spoil the validity of (j23), which is required by gauge invariance. Three- loop 
effects in the Higgs propagators would be taken into account but three-loop effects in 
the couplings T^q+q- would be neglected, resulting in fake large effects. Therefore, 
tree-level Higgs-boson masses have to be used in the propagators for the bosonic 
two- loop contributions (as has been done in Fig. HJ). 

The shift induced by replacing the tree-level Higgs-boson masses by loop corrected 
masses in the two-loop contributions gives an indication of the possible size of bosonic 
three-loop corrections. A shift of O(50 GeV) in irth induces a variation in a M of up to 
about 0.2 x 10~ 10 as shown in Fig. |21 which can be regarded as a three-loop uncertainty. 

5 Chargino/neutralino loop contributions 

The 2-loop contributions to a M containing a closed chargino/neutralino loop constitute 
a separately UV-finite and gauge-independent class. The corresponding diagrams are 
shown in Fig. In this section we discuss the numerical impact and the parameter 
dependence of this class. 

The chargino/neutralino two-loop contributions, a*' 2L , depend on the mass param- 
eters for the charginos and neutralinos fi, M 12 , the CP-odd Higgs mass Ma, and tan (3. 

5 In Ref. [42], eqs. (27,28), a similar relation between the Higgs mass matrix T ! j )i< f )j and the gauge 
dependence of tadpole diagrams has been derived. Gauge dependent tadpole diagrams are in par- 
ticular the ones with G ± loops, thus involving the coupling T^.q+q- , so that this constitutes again 
a gauge relation between T ( j >i< f )j and T^.q+q-. 



13 



\i \i \i 




Figure 5: Generic two-loop diagrams to a M with a closed chargino/neutralino loop. <p, 
ip denote the scalar particles h, H, A , H ± and G 0,± ; V denotes the vector bosons 7, 
Z, W ± ; x stands for the chargino/neutralino mass eigenstates X\2i X? 2 3 4- 



Since we use the loop-corrected values for the CP-even Higgs masses M htH [20,21,27] 
in the propagators for this class of diagrams, formally an effect of 3-loop order, a*' 2L 
also has a slight dependence on other MSSM parameters through Mh,H (we denote 
the loop-corrected masses as M^h and the tree- level masses as ui^h] see also the 
discussion in Sect. HJ). It is interesting to note that, contrary to Ref. [24], no tree- level 
relations in the Higgs sector were needed in order to find a UV- finite result. This is 
due to the fact that each two-loop diagram contributing to (g — 2) M together with its 
corresponding subloop renormalization is finite. 

It turns out that the parameter dependence of a*' 2L is quite straightforward. If 
all supersymmetric mass scales are set equal, 6 /1 = M 2 = Ma = M SU sy, the approx- 
imate leading behaviour of a^ ,2L is simply given by tan/?/M| USY , and we find the 
approximate relation 

v2L 1n ,n 10 /tan /A / 100 GeV\ 2 . , , 

As shown in Fig. the approximation is very good except for very small Msusy and 
small tan/3, where the leading term is suppressed by the small fi, and subleading 
terms begin to dominate. 

The chargino/neutralino sector does not only contribute to a^' 2L but already to 
af 3 SY < 1L , so it is interesting to compare the one- and two-loop contributions. For 

6 For simplicity we assume throughout this paper that M\p are related by the GUT relation 
Ah = 5/3 4/4 M 2 . 
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Figure 6: Comparison of the full result for a*' 2L with the approximation f° r 
\i = M 2 = Ma = M SU sY- The full lines show the complete result, while the dashed 
lines indicate the approximation. For the upper curves tan j3 = 50 has been used, for 
the lower ones tan (3 = 5. 
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Figure 7: Comparison of the supersymmetric one-loop result qSusy.il (dashed) with 
the two-loop chargino/neutralino contributions a^' 2L (dash-dotted) and the sum (full 
line). The parameters are \i = M 2 = Ma = Msusy, tan/3 = 50, and the sfermion 
mass parameters are set to ITeV. 
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the case that all masses, including the smuon and sneutrino masses, are set equal 
to M SU sY, the one-loop and two- loop contributions can be trivially compared using 
eqs. ©, (23), showing that the two- loop contribution shifts the one-loop result by 
about 2%. 

However, the chargino/neutralino sector might very well be significantly lighter 
than the slepton sector of the second generation, in particular in the light of FCNC 
and CP-violating constraints, which are more easily satisfied for heavy 1st and 2nd 
generation sfermions. In Fig. [7| the chargino/neutralino two- loop contributions are 
therefore compared with the supersymmetric one-loop contribution a^ USY,1L at fixed 
high smuon and sneutrino masses Mj = 1 TeV. The other masses are again set equal, 
\i = M 2 = Ma = Msusy- Furthermore, we use a large tan (3 value, tan/5 = 50, which 
enhances the SUSY contributions to a M . 

We find that for Msusy ~ 400 GeV the two-loop contributions become more and 
more important. For Msusy ~ 100 GeV they even amount to 50% of the one-loop 
contributions, which are suppressed by the large smuon and sneutrino masses. 

It is noteworthy that in the numerical analysis of the sfermion loop contributions 
in Ref. [15] it was crucial to take into account the experimental constraints on the 
light Higgs-boson mass [18-21], electroweak precision observables [22-24] and the 
6-decays B — > X s ^ and B s — > ji + fj,~ [25,26]. In the case of chargino/neutralino loop 
contributions, however, these experimental constraints have only little impact on the 
possible choices of //, M 2 and M4 since the relevant observables depend strongly on 
the 3rd generation sfermion mass parameters. 

From Figs. El 13 one can read off that the chargino/neutralino two-loop contribu- 
tions can amount to about ±10 x 10 -10 if //, M 2 , and Ma are around 100 GeV. This 
parameter region is rather constrained but not entirely excluded by the experimental 
bounds from direct searches for charginos and neutralinos, from Higgs searches and 
from 6-physics. 

The dependence of a^ USY,1L + a* ,2L and a^ USY ' 1L on /x and M 2 is shown in the 
contour plots of Fig. |HJ We fix Ma = 200 GeV and vary tan j3 and the common smuon 
and sneutrino mass Mj, which has an impact only on the one-loop contribution. We 
have checked that these parameter choices are allowed essentially in the entire /i-M 2 - 
plane by the current experimental constraints mentioned above, provided the t and 
b mass parameters are of 0(1 TeV). The contours drawn in Fig. |H] correspond to the 
lcr, 2<T, . . . regions around the value a^ xp — a ^ heo > SM = (24.5 ± 9.0) x 10~ 10 , based on 
Refs. [3,7]. We find that for the investigated parameter space the SUSY prediction 
for a M lies mostly in the — 2 a region if \x is positive. However, the new two-loop 
corrections shift the 1 a and 2 a contours considerably. This effect is more pronounced 
for smaller tan /3 and larger Mj. 

The Ma- dependence of the pure two-loop contributions a^ ,2L is shown in Fig. |U] 
for /i, M 2 = 150, 500 GeV. For the smaller values of // and M 2 , a^' 2L behaves ap- 
proximately as I/M4; for the larger values of \x and M 2 the M^-dependence is less 
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Figure 8: Contours of a^ USY ' 1L + a* ,2L (solid border) and a^ USY ' 1L alone (dashed line) 
in the /z-M 2 -plane for M A = 200 GeV. The slepton mass scale (which enters only 
the one-loop prediction) and tan/? are indicated for each plot. The contours are 
at (24.5,15.5,6.5,-2.5,-11.5,-20.5) x 10~ 10 corresponding to the central value of 
a ex P _ a theo,SM = ( 245 ± 9.0) x 10~ 10 and intervals of 1-5(7. 
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Figure 9: Dependence of a*' on the CP-odd Higgs mass Ma for tan/3 = 50 and 
fl = M 2 = 150, 500 GeV. 



6 Conclusions 

We have calculated two kinds of two-loop contributions to (g — 2) M : the purely bosonic 
two- loop corrections (diagrams involving /i, gauge and Higgs bosons) in the SM 
and the MSSM, and the two-loop corrections involving a neutralino/chargino subloop 
in the MSSM. 

In the SM our calculation provides an independent check and a slight exten- 
sion of Refs. [12,13], which was up to now the only calculation of the bosonic two- 
loop corrections, using the approximation M h sm ^> Mw- We find good agree- 
ment with Refs. [12,13]. The variation with M h sm is at the level of 0.3 x 10 -10 
for 100 GeV < M h sm < 500 GeV. The final result can be well approximated by the 
Mfl-sM-independent logarithmic piece, eq. (j22J), a^ os,2L = (—2.2 ± 0.2) x 10~ 10 , where 
the central value corresponds to M h sm m 200 GeV. 

In the MSSM our calculation completes the program begun in Ref. [15] to calculate 
all MSSM two-loop corrections to two-Higgs-doublet model one-loop diagrams. This 
class of corrections has a potentially larger impact on the MSSM prediction of (g — 2) M 
than the higher-order corrections to purely SUSY one-loop diagrams, since it is not 
suppressed in case the smuon and sneutrino masses are large. In addition to the SM 
fermion and scalar fermion loop corrections derived in Ref. [15], we have calculated 



18 



here the diagrams with a neutralino/chargino loop as well as the purely bosonic two- 
loop diagrams in the MSSM (i.e. diagrams involving no SUSY particles). 

For the purely bosonic two-loop corrections in the MSSM we find only a small devi- 
ation from the SM result. The logarithmic piece is identical, while the non-logarithmic 
piece causes only a difference in a M below the level of 0.3 x 10 -10 compared to the SM 
result with a light Higgs. The result for the purely bosonic two-loop corrections in 
the MSSM lies in the range -1.5 x lO" 10 < a £ os < 2L > MSSM < - 2.0 x 10" 10 , depending 
on the values of Ma and tan/3. 

A sizable shift, on the other hand, arises from the two-loop contributions where a 
neutralino/chargino loop is inserted into a two-Higgs-doublet model one-loop diagram. 
Depending on the values of /i, M 2 , M4, and tan/3, the result can amount up to 
±10 x 10~ 10 . For example, if all masses are set equal to a common scale Msusy we 
find the approximate relation 

v2T in /tan /A /100 GeV\ 2 

Thus, the regions in the MSSM parameter space which satisfy the requirement that 
a M agrees with the experimental value at the 2 a level can be shifted significantly as 
a consequence of the new two-loop corrections. 

The analytical results presented here will be implemented into the Fortran code 
FeynHiggs (see: www.feynhiggs.de). Based on the new results presented in this 
paper, the last missing piece of the MSSM two-loop contributions to (g — 2)^ are the 
two-loop corrections to the SUSY one-loop diagrams. Their inclusion, which will be 
desirable in order to further reduce the theoretical uncertainty, will constitute the first 
complete two-loop result within the MSSM. 
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